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ABSTRACT 

Although simple radial equilibrium turbine design 
has long been a useful design method for axial flow gas 
turbine blading, more demanding requirements for high 
performance machinery make such an approximation less 
justified. More exact analytical methods do exist, but 
require many times more effort, even if a computer is 
utilized . 

A relatively simple model for improving a simple 
radial equilibrium preliminary design is proposed in 
this thesis. The axial velocity -distribution, determined 
from the input SRE design, is modified to take into 
account the results of flow studies in flared annuli 
and of actuator disk theory. Using the corrected axial 
velocity distribution, new velocity triangles can be 
calculated. 

A FORTRAN IV computer program to implement the 
proposed model was v^ritten and is presented. Its 
potential is demonstrated by application to three 
example simple radial equilibrium designs. The results 
show practical differences exist between a simple radial 
equilibrium design and the design corrected by the 
proposed method. 

Thesis Supervisor: A. Douglas Carmichael 

Title; Professor of Power Engineering 

Thesis Reader: David Gordon Wilson 

Title: Professor of Mechanical Engineering 
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NOTATION 



A flow area, ft^ 

a coefficient in equations 2.24, 2.25, 2.27, 

and 2.28 

A annulus area, ft^ 

ann 

b coefficient in equations 2.24, 2.25, 2.27, 

and 2.28 

Cp specific heat at constant pressure, BTU/lbin/®R 

c^ specific heat at constant volume, BTU/lbm/°R 

32.2 ft lbra/lbf/sec=^ 

h specific enthalpy, BUT/lbm; actuator disk 

height, ft 

J 778 ft Ibf/BTU 

M Mach number 

m mass flow rate, Ibm/sec 

N turbine rpm, revolutions per minute 

P pressure, Ibf/ft^ 

t 

Q rate of heat transfer to fluid, BTU/sec 

r radial distance from machine axis, ft 

R specific gas constant, ft lbf/lbm/°R 

s specific entropy, BTU/lbm/°R 

T temperature, °R 

u blade velocity, ft/sec 

V fluid velocity, ft/sec 

W rate of output shaft work, BTU/sec 

W relative velocity, ft/sec 

X axial distance from an actuator disk, ft 
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a 

Aa 

B 

AB 

Y 

P 

4 > 

to 



absolute angle measured from axis, degrees 

stator turning angle, degrees 

relative angle measured from axis, degrees 

rotor turning angle, degrees 

specific heat ratio 

density, Ibm/ft^ 

flare angle, degrees 

angular velocity, radians/sec 



Subscripts 



h 

1 

0 
r 

rel 

t 

X 

0 

1 

2 
3 

+ 00 



hub condition 

reference condition 

total (stagnation) condition 

radial component 

relative 

tip condition 

axial component 

tangential component 

rotor inlet; plane one condition 

rotor outlet; plane two condition 

plane three condition 

condition far downstream 

condition far upstream 



NOTE: Notation not appearing here may be computer 

program variables listed in Appendix A. 
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1 . 



INTRODUCTION 



1.1 General 

Energy is a topic of increased interest in current 
times. Searches for new sources of energy and new means 
of converting energy into useful forms have been under- 
taken, along with efforts to conserve the resources 
available. While exploiting solar energy, tapping 
geothermal power, and drilling for oil off-shore receive 
much public interest, progress in improving existing 
conventional forms of energy transfer and conversion must 
also be made. The great importance of energy utilization 
has made it worthv;hile for designers of energy related 
systems to pay even more attention to achieving the best 
possible performance. Improving the efficiencies of the 
various hardware components which carry out the heat 
transfers and energy conversion functions of- the thermal 
cycle is one desirable goal. To accomplish this, it is 
useful to make design methods more accurate and more 
easily applied to real problems. 

1.2 Axial Flow Turbines 

This thesis deals with one commonly used component 
for converting energy into useful work--the axial flow 
turbine. Turbines convert the kinetic energy v/hich has 
been provided the v;orking fluid into mechanical work 
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performed upon a rotating shaft. Axial flow turbines are 
so named because the flow of the \>/orking fluid through 
the machine is predominantly in a direction parallel to 
the machine's axis of rotation. Such turbines can 
operate with a variety of working fluids, and in either 
open or closed thermal cycles. The most common con- 
figurations are as a steam turbine in the Rankine Cycle 
(a closed cycle) or as a gas turbine in the Brayton 
Cycle (an open or closed cycle) . Historically, the 
steam turbine found earlier application than the gas 
turbine, which recently has gained in importance. The 
development of the gas turbine depended upon advances in 
high temperature material properties and was also closely 
tied to the successful operation of axial compressors. 

1.3 Aerodynamic Design of Axial Flow Turbines 
1.3.1 General 

One aspect of the design, which precedes similarly 
for axial turbines of all kinds, is the aerodynamic 
design of the turbine blading. The efficiency of the 
turbine depends upon the smooth transfer of kinetic 
energy from the fluid flow, with as little energy wasted 
as possible. The shapes and angles of the turbine 
blading must be appropriate to the pattern of the flow. 
The aerodynamic design of the turbine determines what 
these flows and shapes should be. 
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In general, the aerodynamic designer of an axial 
turbine follows the steps listed below: 

(1) for an assumed operating point and 
efficiency, carry out the thermodynamic analysis of the 
thermal cycle to determine the mass flow, and inlet and 
outlet states for the turbine; 

(2) postulate a working fluid flow which 
satisfies the laws of motion and produces the required 
work ; 

(3) design the shape of the turbine blading 
appropriately; 

(4) calculate the efficiency and evaluate the 

design ; 

(5) repeat all of the above steps as necessary. 
1.3.2 Flow Solutions 

As is true of most design efforts, there are many 
levels of sophistication and detail which the designer 
may seek, depending upon the stage of the design and/or 
the cost and criticality of the product. In the aero- 
dynamic design of a turbine, the designer must decide 
which level is best in each case. For, although the 
equations of motion and continuity for a fluid are 
readily v;ritten, they cannot be solved in practice unless 
simplifying assumptions are made. Fortunately, such 
simplified forms have proved adequate in many applications. 



13 



1.3.3 Two-Dimensional Flow 



If the hub- tip ratio of a turbine is nearly 1.0, 
no large radial component of fluid velocity is expected, 
and flow conditions will be approximately the same from 
root to tip. The flow can then be analyzed in two 
dimensions. This results in blades of constant cross- 
section, which have been extensively studied experimentally 
in the form of cascades. 

1.3.4 Three-Dimensional Flow 

In some cases the designer faces requirements on 
blade tip speed, pressure ratio, and mass flow which 
could result in low hub-tip ratios (i.e., relatively 
long blades) and considerable turbine flare. These 
conditions x'educe the accuracy of simple two-dimensional 
design methods and force the consideration of three- 
dimensional flov/. The general solution is impractical, 
but simplified methods have been developed and used with 
success. Several of these approaches are mentioned below. 

1.4 Approaches to Three-Dimensional Design 

The first assumption made in virtually all cases is 
to assume axial symmetry. The number of blades is 
sufficiently large that circumferential conditions at a 
station can be represented by average values.^ 

1.4.1 Radial Equilibrium 

Assuming axisymmetric flow, the differential 
equation of motion in the radial direction for the fluid 
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can be put in a form (see paragraph 2.4.1) which may be 
numerically integrated by a computer. The procedure in- 
volves complicated iterations, however, and elaborate 
measures must be taken to ensure the proper convergence 
of the streamline locations. Carmichael has discussed 
the procedure and presented block diagrams of possible 
computer programs.^ This method allows for fluid 
velocities in the radial direction which result from 
imbalances in the static pressure and centrifugal 
acceleration forces on the fluid, but its complexity 
is formidable. 

1.4.2 Simple Radial Equilibrium 

If the radial components of the fluid velocity are 
neglected, the so-called streamline curvature term in the 
differential equation drops out. The differential 
equation of motion in the radial direction can then be 
written involving only the axial and tangential velocities. 
This can be easily integrated for axial velocity when 
given a reasonable distribution for the tangential 
velocity (see. paragraph 2.4.3). Simple pairs of functions 
for axial and tangential velocity, so determined, have 
been widely used in three-dimensional design. 

1.4.3 Actuator Disk Theory 

Simple radial equilibrium analysis depended upon the 
radial fluid velocity being zero before and after each 
blade row. In most cases, changes in the radial distri- 
bution of density result in radial shifts in the 
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streamlines. Therefore, a non-zero radial velocity 
component must exist. 

An alternative to the simple radial equilibrium 
analysis is the actuator disk formulation. Each blade 
row is replaced by a narrow "actuator disk" which provides 
the appropriate sudden change in fluid tangential velocity. 
This device permits an approximate solution to the 
equations of motion, allowing finite radial velocities, 
but still is restricted to annular flow between walls 
parallel to the machine axis. (See paragraph 2.5) 

1.5 Fluid Flow in Flared Annuli 

A frequent design practice is to keep the axial 

velocity approximately constant throughout the machine, 

compensating for changes in fluid density by increasing 

the annulus area. This area change results in "flare"; 

that is, the hub and shroud surfaces are not parallel to 

the machine axis, but are sloped. 

This flare could be expected to cause a general 

change in the axial velocity distribution within the 

annulus, with significant radial velocity components 

2 

appearing near the walls. Carmichael and Pai have 
suggested a relation which can be used to calculate the 
change in axial velocity in flared annuli. 

1.6 Summary 

The preceding discussion was intended to demonstrate 
a need for a relatively simple method to predict the 
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fluid flow distribution in the presence of flare, improving 
upon the accuracy of the simple radial equilibrium model. 
Such a method would hopefully be useable in the preliminary 
design phase with little increased effort, as compared to 
the more complicated streamline curvature approaches. 

A method is proposed in this thesis for "correcting" 
a preliminary simple radial equilibrium turbine design 
for flare and streamline curvature effects. This is 
accomplished by modifying the axial velocity radial 
variation, as determined from the simple radial equilibrium 
equations. A correction for turbine flare is applied 
using the findings of Pai; next, the axial velocity is 
further modified at the blade leading and trailing edges 
by the approximate solutions to the actuator disk analysis. 
As the computations are still iterative in nature, they 
are very suitable for computerization. 

A computer program to modify a simple radial equili- 
brium turbine design, according to the model described 
above, is presented in this thesis. While not theoreti- 
cally deep-seated, it is hoped the methodology may be 
helpful in some turbine design applications. 
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2 . BACKGROUND 



The purpose of this chapter is to provide the back- 
ground for the methodology employed in the computer model 
to be discussed in Chapter 3. The emphasis is placed on 
those assumptions and relations which directly apply to 
the program model, although some additional development 
is included for perspective and to point to possible 
areas of application and further work. 



2.1 Fluid Mechanics 

Axial turbine flow is frequently modeled as incom- 
pressible, at least in the radial direction, so that the 
equations of motion for an incompressible fluid apply. 

It is also known that viscous effects can usually be 
neglected outside the boundary layer. ^ In cylindrical 
coordinates, the equations of motion for an ideal, incom- 

3 

pressible fluid can be written as 
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where V^, Vq , and refer to components of fluid velocity 
in the radial, tangential, and axial directions, respect- 
ively. P is the static pressure. The continuity 
equation is 






9V 

X 

8x 



0 . 



(2.4) 



2.2 Thermodynamics 
2.2.1 Work 

From the first law of thermodynamics, the steady 
flow energy equation, describing the change in state of 
a fluid system in flov;ing through a control volume, can 
be written as"^ 



Q - W = m [ (h + 



V' 



2g J out 
^o 



- (h + 



V 



2g in^ ' 



(2.5) 



« 

where potential terms have been ignored. Q is the rate 
of heat transfer to the fluid, W is the rate of output 
shaft work, ra is the mass flow, V is the fluid velocity, 
and h is the enthalpy. A useful definition is to make 
stagnation enthalpy, h^. 



h 



o 



E h + 



2g J 
^o 



( 2 . 6 ) 



Using equation (2.6), and assuming that the heat transfer, 

• • 

Q, from an axial turbine is small (compared to W) , 

equation (2.5) becomes 
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{ 2 . 1 ) 



-W 






where the 2 and 1 refer to the outlet and inlet states 
respectively. 

Following reference 1, the rate of output shaft work, 
W, is also equal to the product of the angular velocity, 

( 0 , of the blades and the torque caused by the moments of 
the external reaction forces. The torque is given by the 
change in angular momentum of the fluid, or 



Torque = m(r 2 V 



02 






( 2 . 8 ) 



Since 



W = Torque • w 

and the blade velocity, u, is just cor, then 

_W ^ ^ 2^02 ~ ^ 1^81 ^ 2.9 

♦ g J ' 

m ^o 

In axial machines it is a good assximption to make U 2 =Uj^=u 
for preliminary design. 

2.2.2 Gas Lav/s 

The equation of state for a perfect gas can be 
written 
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P = PRT, 



( 2 . 10 ) 



where R is the specific gas constant. The following 
relations between the specific heats are the standard 
ones. 



C - C = R/J 

p V ' 

Cp = f. (2.11) 

2.2.3 The Second Law 

From the second law of thermodynamics, using Gibbs 

Equation and the perfect gas law, the specific entropy, 

u - 1 - 4 - 1.5 

s, can be written as 

Tds = dh - ^ dP. (2.12) 

p J 

2.3 Compressible Flow Relations 

In all but very low speed turbomachines, in which 
velocity, pressure, and temperature changes are very 
small, the effects of compressibility of the working 
fluid become important. Consideration of three-dimen- 
sional compressible flow is extremely complex. Fortunately, 
the changes in fluid density are often not great in the 
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radial direction at a station, and need to be considered 
only in the axial flov7 direction. 

The one-dimensional compressible flow relations for 
an ideal gas appear in many representations and are 

3 

conveniently used. Horlock has a discussion of their 
derivation and summarizes the results. The following 
forms will be useful in discussion later. The meanings 
of the symbols can be found in the notation list. From 
equation (2.6), since, for a perfect gas, h = c^’T, 



T 

o 



T + 






2g JC 
^o p 



(2.13) 



Using (2.11) and dividing by T^, equation (2.13) can be 
rewritten as 




(^) 
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2g R T 
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The Mach number, M, is defined by 



M = 



V 

a 



V 




yRT 



From equations (2.14) and (2.15), 
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(2.16) 
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The pressure-temperature relation for an isentropic 
process can be shown to apply in the compressible case, 
and gives 




From continuity, the flow area. A, is 



A = 



m 

PV 



I 



(2.17) 



(2.18) 



which can be determined from equations (2.14) and (2.17) 
above . 

2.4 Simple Radial Equilibrium 
2.4.1 General 

As indicated in the introduction, when the hub- tip 
ratio of a turbine stage is further and further from the 
value of 1.0 (i.e., relatively longer and longer blades), 
the two-dimensional flov/ design methods become less and 
less adequate. In some marginal cases, three-dimensional 
details of the design are ignored in favor of other 
considerations, for example, maintaining constant blade 
sections. But modern performance requirements ma)ce such 
practices less attractive. The general three-dimensional 
flow equations are difficult to solve, however, and 
simplifying assumptions are frequently made. The most 
general such assumption is to neglect any variation of 
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properties in the circumferential direction. In the 



axial direction, where fluid property changes may be 
large, the one-dimensional compressible flow relations 
are assumed to describe the flow. Finally, in the radial 
direction the equation of motion in a simplified form-- 
termed the simple radial equilibrium (SRE) equation — is 
often applied. A discussion of the development and use 
of the simple radial equilibrium equation follows. 

2.4.2 The Simple Radial Equilibrium Equation 

Assuming axial symmetry, the equation of motion in 
the radial direction, equation (2.1), can be written as 
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From the identity and by differen- 

tiating with respect to r, equation (2.6) becomes 
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( 2 . 20 ) 



and equation (2.13) can be rewritten in the form 



1 3P _ 3s 3h 

pJ 3r 3r 3r‘ 



Substitution of equations (2.20) and (2.21) into equation 
(2.19) gives 
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At this point 3s/9r and 9h^/9r are often assumed to be 

zero. According to Horlock^ these two assumptions, 

which prove to be nearly true in practice, hold for the 

flow between blade rovv-s of a reversible turbine in which 

delivered work and total pressure drop across the blade 

are the same at all radii. The term, V • 9V /9x, is 

X r 

called the streamline curvature term and is the cause of 
the difficulty, mentioned in the introduction, in 
obtaining a solution to equation (2.22).^ In addition 
to assuming 3s/9r and 9h^/9r are equal to zero, if it is 
assumed that equals zero at the leading and trailing 
edges of a blade row, the streamline curvature term drops 
out and the partial derivatives can be replaced by 
ordinary derivatives, giving 



V 

X 



dV 

X 

dr 



V 

'^0 dr r ' 



(2.23) 



Equation (2.23) is called the simple radial equilibrium 
equation. 

A more direct derivation of the simple radial 
equilibrium equation, which may be more intuitive, can be 
developed from consideration of an infinitesimal volume 
of fluid in equilibrium with the external forces acting 
on it.^ The details are shown in Figure 2-1, on the 
next page. The substitution for the g^/p(dP/dr) term 
precedes as above. 
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Sum of forces in the radial direction = 0 



r-de(P 



dr 



dr) - P-r-d©- — (pr-d9'dr) = 0 



dr 



v| 

r-de-dr 



r • d© • dr 



i ^ ^ 

p dr r 



Figure 2-1 
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In summary, the simple radial equilibrium equation 
is valid for the following conditions: 

(1) ideal, incompressible fluid flow 

(2) axial symmetry 

(3) perfect gas law holds 

(4) 0 before and after blade rows 

(5) dh^/dr = 0 

(6) ds/dr = 0 

2.4.3 Solutions to the Simple Radial Equilibrium Equation 
From equation (2.9) it can be seen that the work 
performed by an axial turbine is related to the tangen- 
tial velocities before and after the rotor, Vq, and Vq„, 

01 V z 

respectively. A design problem for an axial turbine is 
to propose a reasonable tangential velocity distribution 
which will give the desired v;ork, and then by means of 
the simple radial equilibrium equation find the axial 
velocity, V , so that the blade shapes can be determined. 
For simple functions of Vg = Vg (r) , equation (2.23) is 
easily integrated directly for V = V^(r). 

It would appear that independent functions for Vg^ 
and Vg 2 could be chosen and corresponding functions 
V^j^(r) and ^^.2 calculated from the simple radial 
equilibrium equation. However, the Vg ' s are related by 
the work relation, equation (2.9), and by the assumption, 
used in deriving the simple radial equilibrium equation, 
that the stagnation enthalpy, h^, after each blade row, 
does not vary with radius. 
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A familiar, useful set of tangential velocities 
given by Carmichael and Lewis, and described by Horlock^, 
are the following; 
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(2.24) 
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(2.25) 



where a, n, and b (b > 0) may be chosen by the designers. 
Then, from equation (2.9) 
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and using equations (2.24) and (2.25) 



W _ 2TTrN ,2b. _ 47TNb 

60g J ^ r ^ 60g J* ^ 

m -’o ^o 

Equation (2.26) shows the work done (hence h^) is 
independent of radius and has the proper sign for a 
turbine . 

A modification of the set of tangential velocities 
(2.24) and (2.25) is considered in this thesis. Suppose 
a set of tangential velocities of the following form is 
given : 
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(2.27) 
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2 2 
a„r + a,r +a + — 
2 lor 



(2.28) 



Here also the work is independent of radius, for 
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(2.29) 



Here must be greater than if work is to be done by 
the turbine. The constant b^^ is virtually always positive, 
and is usually negative or zero. The advantage of 
allov/ing b^^ and to be different is to be able to 
provide for special cases such as -free vortex, axial 
outflow (a^, a^, a^, b^ = 0 ; b^ 7 ^ 0 ) . 

Substituting an expression of the form of (2.27) 
into the simple radial equilibrium equation and integrating 
from reference point V . = V (r.) to point V = V (r) , the 
following expression for axial velocity is obtained (the 
details are in Appendix D) : 



V 2 = V .=^-1 a,Mr‘*-r.‘*) - ^ a,a,(r^-r ') 



xi 2 2 



3 12 



(4a a„ + 2a. ^ ) (r^-r . ^ ) -6 (a a +a_b^ ) (r-r . ) 
OZ 1 1 Olzl 1 



+ + 2a^b^(p - |-) 

i i 



(2.30) 
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Equations (2.27) and (2.30) for tangential velocity and 
axial velocity form a set which satisfy the simple radial 
equilibrium equation (2.23). It is important to note 
that many factors and various criteria for good practice 
influence the detailed choice of the constants in the 
expressions for Vg . 

Tv/o well knov/n special cases of equations (2.27) 
and (2.28) are the constant reaction Vg distribution and 
the exponential Vg distribution. These tangential 
velocity distributions are presented in Appendix D 
together with the corresponding function for axial 
velocity. 

2.5 Actuator Disk Theory 
2.5.1 General 

It is recalled that the basis for the simple radial 
equilibrium design method is that the radial component of 
the fluid velocity be zero at the leading and trailing 
edges of a blade row. This assumption results in stream- 
lines of the form illustrated in Figure 2-2 (a). A more 
general method of analysis, which takes some account of 
radial shifts in the streamlines between blades, is the 

g 

actuator disk model. As described by Horlock , a rotor 
or stator rov; is replaced, conceptually, at approximately 
midchord, by a thin disk across which the tangential 
velocity change takes place. Such a device allows the 
equations of motion, including radial velocity terms, to 
be solved analytically, although the formulation is still 
complex. 
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2.5.2 Approximate Solution 

With some simplifications, an approximation to the 
complete actuator disk results can be obtained. Radial 
velocities are assumed to be zero at the walls and small 
elsewhere, and simple radial equilibrium is assumed to 
be valid a sufficient distance from the disk. According 
to this model, the streamlines v^;ill be similar to those 
shown in Figure 2-2 (b) . The approximate solutions for 
the axial velocity near a single disk, in terms of the 
axial velocity on the same streamline far up and down- 
stream, are: 



V (X) =(V ) 

X X -<» 



(V ) - (V ) 

. X +<» X ,'TTX. 

+ exp , 



(2.31) 



upstream of the disk (x < 0) , and 



V (x) = (V ) , 

X X +“ 



(V ) - (V ) 

exp — yp) , 



(2.32) 



downstream of the disk (x > 0). The notation (V ) 
and (V^.)_^ refer to the axial velocities given by simple 
radial equilibrium a long way downstream and upstream, 
respectively; h is the height of the actuator disk. 

This analysis, like the simple radial equilibrium 
analysis described in the previous section, treats ideal, 
incompressible flow in cylindrical annuli. 
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For real turbines, where many rows of blades (disks) 
may influence the axial velocity distribution, the effects 
of each disk may be considered separately, then super- 
imposed to find the real axial velocity.^ This procedure 
is used in Appendix E to find the axial velocity distri- 
bution at a point influenced by two nearest adjacent 
disks . 

2.6 Flow in Flared Annuli 
2.6.1 General 

In addition to improving the methods for cylindrical 
flow, efforts to understand the fluid flow in ducts of 
changing annular cross-section have been made. Such 
sections are commonly used in practice. 

Korlock, in reference 6, describes some theoretical 
models, proposed by various workers, for the determination 
of axial velocity. Analyses by Lewis and Horlock for 
incompressible flow in conical ducts, by Wu for compressible 
flov; in flared annuli, and a simplification of Wu's analysis 
by Walker are discussed. The drawback of these analytical 
methods is their complexity. 

2 

Experimental work by Carmichael and Par has 
suggested a simple expression relating the axial velocity 
of flow in a flared duct to the angle the streamline 
makes with the axis. Their work, described below, forms 
the basis of the flare correction to axial velocity 
given in this thesis. 
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2.6.2 Experiments Representing Fluid Flow in Flared 

2 

Turbines , 

Experiments were conducted using a tilted-floor 
electrolytic tank to measure flows within a flared 
annulus. An electrolytic tank provides an analogue to 
ideal incompressible fluid flow. Electrical potential 
measured in the tank is analogous to fluid potential and 
the electrical potential gradient is analogous to the 
fluid velocity. A sector of an axisymmetric model can 
be represented by tilting the floor of the tank. 

Although flow in a real turbine is compressible, 
the radial variation in density at any axial station is 
relatively small. Consequently, it can be assumed that 
in usual geometries the streamline positions for the 
real, compressible fluid are vei*y nearly the same as for 
the ideal, incompressible case. 

One aspect of the fluid flow measured was the 
general change in axial velocity distributions in the 
flared annulus. Three models were studied, v/ith 
different distributions of flare between the hub and 
shroud walls; each model had a hub-tip ratio of 0.9. 

The distribution was expected to be a function of location 
within the annulus, described by an angle, (p , defined as 
shown in Figure 2-3. The flow conditions were measured 
at an axial position, free of discontinuities, such as 
section A-A' in Figure 2-4. The axial velocities were 
found to have the distribution shown in Figure 2-5, 



34 




Figure 2-3 





axis 

Figure 2-It 




Figure 2-5 



36 



approaching the flow from a source. The equation, 
according to potential flow theory, for the axial 
velocity distribution for a two-dimensional source is 



V 

/ - ((j)=0) = (2.33) 

Pai suggests that comparison of this equation, also 
plotted in Figure 2-5, with the experimental results, 
indicates that for angles of (f> < ± 15°, equation (2.33) 
is a good representation of the axial flow. Even at 
larger angles, the agreement may be good enough to use 
equation (2.33) for some design applications. 

2.7 Axial Turbine Flow Geometry 
2.7.1 Velocity Triangles 

It should be restated that the purpose of deter- 
mining the details of the fluid flow is to be able to 
design the blades appropriately. The relevant fluid 
velocities and flow angles are conveniently displayed 
in the form of vector diagrams, often called velocity 
triangles. Velocity triangles for a single stage axial 
turbine with pure axial inflow are shown in Figure 2-6, 
v;here any radial velocity components have been neglected. 
The calculation of such velocity triangles for each 
leading and trailing edge and at several radii must be 
accomplished before the detailed design of the turbine 
blading can continue. 
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2.7.2 Annulus Area 



Suppose a fluid is flowing uniformly through a small 
area at a given rate. If the area encompassing the flow 
is perpendicular to the streamlines, it is termed the 
flov; area. A, and is given by equation (2.18). If the 
area is in a surface whose perpendicular makes an angle 
a to the streamlines, as in Figure 2-7, then the area 
required to encompass the same flow is 

A = A/cos a. (2.34) 

3.nn 

In the case of an axial flow turbine, A represents 

a.nn 

the annulus flow area perpendicular to the machine axis, 
where a is the same a as in Figure 2-6. From Figure 2-7 
it can be seen that 



cos a = V /V; 

v/hence equation (2.18) becomes 



A — ..TT • 

ann pV 



(2.35) 
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Figure 2-»7 
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3. MODEL DESCRIPTION AND IMPLEMENTATION 



3.1 General 

In the introductory chapter the need for a simple 
methodology to generate a preliminary axial turbine design, 
taking account of turbine flare and radial streamline 
shifts, v;as discussed. The simple radial equilibrium 
design method was presented in Chapter Two, together 
with the results of actuator disk theory and axial flow 
distribution studies in flared annuli. In this thesis 
a procedure is proposed which accepts a simple radial 
equilibrium preliminary design as a basis, and then 
modifies the axial velocity distribution in accordance 
with the actuator disk and flare corrections mentioned 
above. These modifications to the simple radial equili- 
brium flov; solution result in slightly different blade 
shapes, which should be more commensurate with the real 
flow. 

This chapter will present the model for modifying a 
simple radial equilibrium preliminary design, and describe 
a computer program v;hich implements the model. 

3.2 Model Description 

To begin, it is useful to consider a sketch of the 
geometry and then a bare outline of the problem. 
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Figure 3-1 represents a single stage axial flow turbine. 
The stator inlet, rotor inlet, and rotor outlet planes 
are labeled 1, 2, and 3 respectively. Other features 
will be discussed as they become important. 

The following^ are steps in the procedure which may 
be identified in order to guide the discussion: 

(1) "input" a simple radial equilibrium 
preliminary design 

(2) at each of the three principal planes of the 
turbine stage, calculate a set of streamlines which is 
consistent v.’ith certain given conditions and the 
modified axial velocity 

(3) at every streamline in each principal 
plane, calculate velocity triangle data. 

The rest of this section will discuss the above steps in 
general terms. The details of the procedure are best 
understood in the context of the computer program des- 
cription, beginning in paragraph 3.3. 

3.2.1 The Input Design 

The purpose of the procedure discussed in this 
thesis is to make slight "corrections" in an existing 
simple radial equilibrium design in order to refine the 
detailed blade shapes. As a result of carrying out a 
preliminary design for a single stage, many parameters 
will have been established. The working fluid, total 
temperatures and pressures, hub and tip radii, machine 
rpm, mass flov;, and blade widths are known; also, in the 
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case of simple radial equilibrium, an appropriate pair of 
functions for axial and tangential velocity at each plane 
will have been picked. These quantities are indicated in 
Figure 3-1. 

For present purposes, it is assumed that it is 
desired to retain the input annular geometry (hub and 
tip radii, and blade widths) in the modified case, 
although other choices are possible. For example, the 
average axial velocity could be fixed, and the tip radius 
allowed to change to accommodate the flow. 

For simplicity, the input preliminary design is 
restricted to having hub and shroud surfaces of uniform 
slope over the v/hole stage, as in Figure 3-1. Also, the 
rotor tip clearance and stator-rotor gap are assumed to 
be zero. The model could be applied to more general 
cases v.’ith relatively minor modifications to the computer 
program. The turbine flare may be distributed between 
the hub surface and shroud surface in any proportion. 

It should be noted that the procedure described 
here will, in no sense, improve a bad preliminary design. 
3.2.2 Calculating the Streamlines 

To fully determine the blade shape, fluid flow 
information, including axial velocity, must be known at 
several radii. If the design were merely simple radial 
equilibrium, the velocity triangles could be independently 
found at any convenient radius, since and Vg (and U) 
are knovm functions of r alone. However, in the present 
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case, the actuator disk correction to the axial velocity 
at a location on one plane is dependent upon the axial 
velocity along the same streamline at stations up and 
downstream. For this reason a set of streamlines passing 
all three principal planes must be calculated. 

To find the final set of streamline positions, three 
different iterative loops are involved. They are: 

(1) at a given plane and streamline, iteration 
is required to obtain the next (radially outward) 
streamline position; 

(2) at a given plane, iteration is required to 

obtain an average axial velocity ivhich v;ill result in a 

match between the fixed mass flow rate and annulus area; 

(3) since a change in the axial velocity at a 

streamline in one plane affects the axial velocity in 

another plane, the streamlines are recalculated at each 
plane in succession until the changes from one iteration 
to the next become insignificant. 

The modified axial velocity, used in the streamline 
calculations, is found utilizing the approximate actuator 
disk results described in paragraph 2.5 and Appendix E, 
and the flared annulus axial flow relation described in 
paragraph 2.6. The only changes to the simple radial 
equilibrium preliminary design are a result of these two 
"corrections" to the axial velocity, which impact on the 
velocity triangles. The changes represent a departure 
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from simple radial equilibrium. The justification for 
this procedure is that the adjustments are small and are 
improvements in the direction of the real solution. 

3.2.3 Velocity Triangle Data 

After the positions of streamlines have converged, 
all relevant parameters in the velocity triangles can be 
calculated from the primary variables U, corrected V , 
and Vq, and input parameters. Other useful design data, 
for example, stator and rotor fluid turning angles and 
W 2 /W 1 ratio, are easily derived from the basic velocity 
triangle data. 

For purposes of comparison, sets of velocity tri- 
angles can be computed at the same streamline-plane 
locations, but using the unmodified simple radial 
equilibrium axial velocity. 

3.3 Computer Implementation of the Model 
3.3.1 General 

The iterative nature of the detailed calculations 
virtually precludes hand calculations of the streamlines 
as a practical matter. A computer program to accomplish 
the calculations has been written in FORTRAN IV language 
making use of the model easy and inexpensive. 

The goal of the work described in this thesis was 
to develop the model and to program it, so that it 
could be applied in future design studies. Since the 
model itself is experimental, the program is not merely 
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the mechanization of standard, accepted practice. For 
this reason, it was thought advisable to include here a 
sufficiently detailed description of the computer program 
so that the model used is clear and so that changes can 
be made easily in the future if desired. References 7 
and 8 were useful as programming language guides. 

The programming philosophy was to simplify the 
understanding and use of the program and to provide 
flexibility. The program was segmented into subprograms, 
each having a specific function. Identical variable 
labels have the same meaning in all segments of the 
program in order to simplify understanding it. The use 
of COMjMON reduces the storage demands of the program. 
Extensive execution information, such as the number of 
loops performed and intermediate values can be printed 
out at the user's option. Diverging iterative loops are 
terminated by internal testing. 

The following paragraph lists the significant steps 
in the execution of the program to modify a simple radial 
equilibrium design in accordance with the model. Sub- 
sequent paragraphs will discuss each step in detail. 

3.3.2 Procedural Steps 

The following list contains a list of all significant 
steps in the program. Some details, such as initializa- 
tion of variables and testing for loop convergence, are 
omitted and are self-explanatory or can be conveniently 
discussed separately. Free use will be made of the 
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notation employed in the computer program itself. The 
list is intended to serve as a guide for the discussion; 
full appreciation of the interrelationships can be gained 
only by studying the program flow charts. 

(1) Read and process input data 

(2) At the first plane calculate a set of 
streamlines satisfying continuity using 
the simple radial equilibriim axial 
velocity corrected for flare 

a) the first streamline follows the hub 

b) guess at the location of the next 
streamline and iterate until continuity 
is satisfied between two streamlines 

c) find locations of succeeding streamlines 
in turn until the last one has been 
calculated 

d) compare the radial location of the last 
streamline to the desired tip diameter 
at the plane under consideration 

e) if a discrepancy exists, modify the 
"average" axial velocity at that plane 
and repeat from a) above 

(3) Repeat (2) for planes 2 and 3 

(4) After a set of streamlines is calculated 
for all three planes, recalculate them all 
as before but using the actuator disk 
correction on the axial velocity as v/ell 
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(5) Test to see if the streamline positions 
for plane three have changed much; if they 
have, go back to step (4) and recalculate 
the streamline positions for all three 
planes 

(6) Calculate and print velocity triangle data 
for each streamline position at each plane, 
using the modified axial velocity 

(7) Repeat (6), but using the simple radial 
equilibrium axial velocity. 

3.3.3 Overall Program 

The overall program is represented by the flow 
chart in Figure 3-2, on the next page. The circled 
numbers provide points of reference for more specific 
flow charts in later figures. The initialization, 
incrementing, and printing of Ll, L2 , and L3 are included 
in the flow charts so that the user can readily interpret 
certain diagnostic output to be described later. 

3.3.4 Reading and Processing Input Data 

The first step in the program is to read a set of 
input data and perform some preliminary calculations. 

This is accomplished by the main program, f-lAIN. 

The. following is a list of the input, determined 
from a simple radial equilibrium design, which must be 
submitted: 

(1) working fluid, specific gas constant, and 
specific heat ratio 
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(2) 


turbine mass flow rate and rpm 


(3) 


total temperatures and pressures at all 
three planes 


(4) 


hub and tip radii at planes one and three 


(5) 


stator and rotor widths 


(6) 


coefficients of terms in the simple radial 
equilibrium relation for tangential velocity 
at each plane 


and ( 7 ) 


a reference velocity and radius used to 
find the simple radial equilibrium axial 
velocity at each plane. 



In addition, the mimber of streamlines desired and 
certain test values are read in. 

The follov/ing calculations of a preliminary nature 
are performed: 



(1) 


test to ensure the overall flare is not 
zero and calculate flare angles 


(2) 


calculate the basic parameters used to 
apply the flare correction to the axial 
■velocity (paragraph 3.3.9) 


(3) 


calculate dimensional test values from 
input (paragraph 3.3.13) 


(4) 


calculate blade heights 



(5) calculate certain other often used combina- 
tions of input which remain constant for 
the data set under consideration 
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Selected parameters are printed out prior to further 
execution of the program. 

3.3.5 Calculating a Streamline Location 

The heart of the whole program is the calculation of 
streamlines; the most basic step in the streamline 
generation is the calculation of the "next" streamline 
location, RNEXT, from the position of an adjacent stream- 
line, R(I,J), and values VXZ and VTZ, representing average 
values of the axial and tangential velocities in the 
channel between RNEXT and R(I,J). This calculation is 
performed by subprogram RNEXT(I), whose flow chart appears 
in Figure 3-3. 

The total velocity, V, is determined from the axial 
and tangential components, and with equation (2.13) 
determines the static to total temperature ratio. The 
temperature ratio and the equation of state for a perfect 
gas, equation (2.10), are used in another compressible 
flow relation, equation (2.17), to find a density RHO. 

Then the continuity equation (2.35) can be used to find 
a flow area, AREA. From 

RNEXT = +R(I,J)2 (3.1) 

/ 7T 

the location of the next radially outward streamline can 
be found . 
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Figure 3-3 



since RNEXT depends upon average values of axial 
and tangential velocities, which are in turn functions 
of radius, iteration is required. This iteration is 
accomplished by the subroutine STREAM shown in Figure 3-4. 

STREAM utilizes RNEXT as a guess for the next 
streamline and then calculates the average velocity 
values. The average values VXZ and VTZ are then used 
to calculate a new RNEXT and so on. When finally the 
difference between any RNEXT and the previous RNEXT is 
small, STREAM returns the values R(I,J+1) = RNEXT, 
VX(I,J+1), and VT(I,J+1) to the main program. The index 
L3 counts the number of times RNEXT is calculated during 
the calculation of a set of JJ streamlines at any plane. 

In addition to finding an R(I,J+1), STREAM also 
applies the flare and actuator disk corrections to the 
unmodified (simple radial equilibrium) axial velocity. 

The corrections are accomplished by multiplying the 
unmodified axial velocity, calculated in subprogram FVX(I) 
(paragraph 3.3.6), by the corrections for flare and 
actuator disk .effects which are calculated by subprograms 
FLARE(I) and DISK(I), respectively. (See paragraphs 
3.3.9 and 3.3.10) Since the very first time a set of 
streamlines is calculated at a plane, DISK is undefined, 
only the flare correction can be utilized in the first 
(Ll=l) loop. 
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3.3.6 Calculating and Vg (Simple Radial Equilibrium) 
In a simple radial equilibrium preliminary design, 

it is standard procedure to pick a functional form for 
the radial variation of tangetial velocity. Then, the 
radial variation of the axial velocity is determined by 
the simple radial equilibrium equation (2.23). In the 
present case, the tangential velocity can have the form 
of equation (2.27) where the coefficients a^, a^, a^, and 
b are chosen by the designer. The simple radial equili- 
briiom equation has been integrated for this form in 
Appendix D, resulting in equation (2.30). The function 
subprograms FVX(I) and FVT(I) calculate the sim.ple radial 
equilibrium values of axial velocity and tangential 
velocity, respectively, given an argument RR and input 
parameters A2(I), A1(I), A0(I), B(I), VX0(I), and R0(I). 

3.3.7 Generating a Feasible Streamline Set at a Plane 
The main program calculates a set of streamline 

locations at a plane starting with the hub and using 
STREAM to step outward to the tip. Because the axial 
velocity has been modified, however, the average axial 
velocity is different from the simple radial equilibrium 
design value, and the last streamline does not fall on 
the fixed tip radius. To correct this, a reference 
axial velocity, VXR(I), used by FVX(I), is adjusted and 
the streamline set is recalculated. 
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The iteration is continued until the last streamline 
location is close enough to the desired tip radius. A 
flov; chart of the section of MAIN that generates a set of 
streamlines that fit the geometry is shown in Figure 3-5. 
The variables S and SI do not play a direct role in the 
streamline generation and will be explained in paragraph 
(3.3.11). The adjustment of the reference axial velocity 
is explained in the next paragraph. The index L2 counts 
the number of iterations made at a plane to converge to 
the proper R(I,JJ). At the users option (by specifying 
KK=1 in the input data) L2 and other parameters may be 
printed out after each trial set of streamlines. 

3.3.8 Adjusting VXR(I) 

Initially, VXR(I) is set to VX0(I), the value 
determined from the simple radial equilibrium design. 
Whenever the location of the last streamline, R(I,JJ), 
misses the tip radius, RT(I), VXR(I) is adjusted to 
increase or decrease the average axial velocity before 
recalculating a new set of streamline locations, in order 
to make R(I,JJ) converge to RT(I). 

To find an expression for a correction to the 
reference axial velocity in terms of the error, 

DR = R(I,_JJ) - RT(I), an expression for the flow 

51 = V -A (3.2) 

P X 
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can be differentiated, giving 

d(^) = AdV^ + V^dA. (3.3) 

But we desire d(m/p) = 0 and A = ir(RT(I)^ - R(I,1)^) 
whence 

2V RT (I ) • DR 

dV = — . (3.4) 

RT(I) ^-R(I,1) 2 

From experience, to make the convergence more rapid, the 
following term was added to the VXR(I) by MAIN: 



dV 



X 



2 VX(I,JJ) DR 
RT(I) - R(I,1) 



(3.5) 



3.3.9 Correcting Axial Velocity for Flare 

The simple radial equilibrium axial velocity 
calculated by subprogram FVX(I) is corrected by multiplying 
it by a non-dimensional factor computed by subprogram 
FLARE(I). The basic relation involved is equation (2.33), 
rewritten as 



FLARE = cos^({) 



1 

/ 

tan^<}) + 1 



(3.6) 



where 



tan ({> 



RR - REF 
D 



(3.7) 
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RR is the radial location of the point of interest, D is 
a parameter depending upon the plane under consideration, 
and REF is a geometric parameter. For additional details 
see Appendix F . FLARE is always positive and less than 
or equal to one. 

3.3.10 Correcting Axial Velocity for Actuator Disk 
Effects 

The actuator disk correction is applied after the 
flare correction and is also a multiplicative non-dimen- 
sional factor, calculated by subprogram DISK(I). The 
relations used in DISK(I) are forms of equations (E.l), 
(E.5), and (E.2) for planes one, two, and three, 
respectively. DISK is always positive and may be 
greater than one. 

3.3.11 Testing for Overall Streamline Convergence 

The corrected axial velocity at a streamline depends 
upon its location and the value of axial velocity on the 
same streamline at planes up and downstream. Therefore, 
in general it is necessary to recalculate the streamlines 
at each plane several times, the streamlines moving less 
with each iteration. The streamline positions are 
calculated in order at planes one, tv;o, and three at 
least twice. Each time the streamlines for plane three 
are recalculated, the following sum is accumulated: 

JJ 

^ ~ T -1 ^^^^'‘^^this iteration ^^^'*^^last iteration 

(3.8) 
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Refer to Figures 3-2 and 3-5. If S is small enough, the 
streamline positions have not changed much from the 
previous iteration. This means that the set of corrected 
axial velocities, VX(I,J), and streamline positions, 
R(I,J), have converged to mutually consistant values, 
and iteration can be terminated. 

3.3.12 Velocity Triangle Calculations 

Once the streamlines have been located at each 
principal plane, a velocity triangle can be constructed 
at each point. PI, P2, and P3 in Figure 3-1 are the 
three points for one particular streamline. JJ triangles 
are constructed at each plane. 

The fundamental input variables to the velocity 
triangle calculations are the streamline location R(I,J), 
VT(I,J) , and VX(I,J), which were calculated concurrently 
in the iteration process. The following derived 
quantities were calculated by MAIN at each plane: 



U = 2-n • RPM • R(I,J)/60 



(3.9) 



V = /VT ( I , J ) 2 + VX(I,J)2 



(3.10) 



W = /VX(I, J) 2 + [U - VT(I, J) ] 2 



(3.11) 



a 



arctan ( 



VT(I,J) 

VX(I,J) 



(3.12) 



B = arctan ( 



VT(I,J) - U 
VX(I, J) 



(3.13) 
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M = from equation (2.16) 



M 



rel 




(3.14) 



In addition, the following quantities were calculated for 
each streamline: 



VI ^ V (plane 2) 
V2 V (plane 1 ) 



W2 _ W (plane 3 ) 
V71 “ W (plane 2) 



stator turning angle = a (plane 2) - a (plane 1) (3.17) 



rotor turning angle = 3 (plane 3) - 3 (plane 2) (3.18) 

To compare the above results to those attainable by 
simple radial equilibrium, all of the above velocity 
triangle calculations are repeated at the same locations, 
but with the unmodified simple radial equilibrium axial 
velocity. 

Figure 3-6 is a flow chart describing the calculation 
and printing of the above results. 

3.3.13 Iteration Test Values 

For flexibility, the user supplies the test values 
used to terminate the three iteration loops. The three 
loops are (1) overall streamline convergence. 
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I^IN: Calculate and Print Velocity Triangles 




Figure 3*6 
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(2) convergence of the J J ' th streamline at a plane to 
RT(I), and (3) the location of the next streamline 
R(I,J + 1). These three loops are counted and labeled 
by the variables LI, L2 , and L3, respectively. Test 
values TESTl, TEST2 , and TEST3 are calculated from input 
param.eters Tl, T2, and T3 and the following relations: 

TESTl = Tl (RT(2) - R(2,l)) (3.19) 

TEST2 = T2 (RT(2) - R(2,l)) (3.20) 

TEST3 = T3 (RT(2) - R(2,1))/JM1 (3.21) 

An iteration loop is terminated normally when the error 
term becomes less than the appropriate test value. 

3.4 Convergence Checks 

If for some reason the program ever begins to 
calculate diverging values of the error term in any of 
the three iteration loops, the program will automatically 
terminate the calculations on the current data set. This 
is accomplished, as depicted in Figure 3-7, by exiting if 
the error term has increased at least twice in that loop. 

3.5 Program Verification 

In this thesis a model to improve simple radial 
equilibrium design by applying flare and actuator disk 
corrections to axial velocity is described and a 
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Test for Loop Divergence 




Fipure 3-7 
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computer program to implement the model is presented. 

The computer program has been verified to the extent that 
there is confidence that it accurately reflects the model. 
Full evaluation of the model, including comparison to 
more complicated, exact design methods, must be left for 
future efforts. 

Although the iterative nature of the streamline 
calculations makes hand calculation of the whole procedure 
very burdensome, it is easy to verify a solution by hand. 
The outputs of subprograms FLARE(I), DISK(I), FVX(I), 
and FVT(I) were all verified by hand, given a computer 
calculated streamline position R(I,J + 1). With these 
values and R(I,J), RNEXT(I) and STREAM Vv’ere verified by 
hand calculating R(I,J + 1). The ’computer calculations 
of flare geometry, velocity triangle related output, and 
other algebraic relationships, all agreed with hand 
calculation . 

Since each input data set must be a complete simple 
radial equilibrium preliminary design, the number of 
independent geometries, working fluids, tangential 
velocity distributions, etc., tested was not very large. 
Nevertheless, it is felt the verification of the program 
was sufficient to ensure accuracy for feasible designs. 
Results are shown for three sample cases in Chapter 4. 
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4 . RESULTS 



4.1 General 

This chapter presents the results of applying the 
computer program to three sample simple radial equilibrium 
designs. The computer output listing is included and 
plots of representative variables are provided for easy 
comparison between the simple radial equilibrium input 
design and the modified design. 

4.2 Free Vortex Design 

Table 4-1 shows the characteristics of an axial 

inflov; and outflow free vortex stage designed as the last 

stage of a double flow turbine in a geothermal power 
9 

application. A sketch of the geometry appears in 
Figure 4-2. Notice the design has a constant hub radius, 
but has large flare in the shroud. The working fluid is 
Freon 21. Eleven streamlines are called for. 

The computer printed results are listed in Tables 
4-3, 4-4, and 4-5. It is seen that the large changes in 
axial velocity (about 25% at the hub) can result in a six 
degree change in the relative inflow angle to the rotor. 
Less pronounced, but readily noticed changes are shown 
in other flow angles, velocity ratios, and Mach numbers. 
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In Figure 4-6, the axial velocity, stator outlet 
angle, rotor fluid turning angle, and relative velocity 
ratio are plotted for both modified and original simple 
radial equilibrium designs. It may be significant to 
note that the changes at the hub (where <p = 0°) may be 
considered to be "indirect" in that they result from the 
increased average axial velocity required to obtain the 
same mass flow. 

4.3 Constant Reaction Design 

The input parameters for another example having the 
same annular geometry as above but with a different 
assumed tangential velocity distribution are shown in 
Table 4-7. Axial inflow to the stage is maintained, 
but the constant reaction type simple equilibrium 
relations (see Appendix D) are used at the rotor leading 
and trailing edges. 

The computer printed results are shown in Tables 
4-8, 4-9, and 4-10, and plots of axial velocity and rotor 
relative flow angles appear in Figure 4-11. From Figure 
4-11, the modification to axial velocity doesn't seem as 
significant as in the first example; the influence on the 
rotor flov; angles is still significant, however 
(A3 = 6° at the hub) . 
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4.4 Free Vortex, High Mach Number Design 



To demonstrate the adequacy of the computer program 
to calculate in regimes where M > 1.0, the preliminary 
design summarized by Table 4-12 was run. In this case, 
the flare is symmetrically distributed between hub and 
shroud, as is shown in Figure 4-13. 

The detailed results are listed in Tables 4-14, 

4-15, and 4-16. The axial velocity and Mach number at 
the rotor inlet (plane 2) are plotted in Figure 4-17 
for comparison with the simple radial equilibrium values. 
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5 . CONCLUSIONS 



5.1 General 

A model for improving a simple radial equilibrium 
design by modifying the axial velocity has been discussed. 
The computer program presented here provides a means to 
quickly and inexpensively incorporate the model into 
preliminary design procedure. 

The improved preliminary design model is basically 
a perturbation of the simple radial equilibrium design 
procedure, and comparisons of the two procedures can be 
made easily from the computer program output. It v/ould 
be useful to compare the improved model results to designs 
produced by more rigorous and complicated analytical 
procedures as v;ell, to see where the balance of complexity 
versus accuracy lies. Since the program described here 
is only a "design tool" and was not derived analytically, 
comparisons to more exact procedures are really needed to 
give full respectability to the results. 

The changes made as a result of applying the model 
to a simple radial equilibrium preliminary design are of 
a magnitude which suggests that they may have measurable 
effects upon turbine performance. Further study is 
needed to quantify these effects in terms of efficiencies. 
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5.2 Improving the Model 



Aside from possible improvements in the generation 
of the correction to the simple radial equilibrium axial 
velocity, other changes in the program to extend its 
usefulness may be desirable. An example is extending the 
range of the program to include more than one stage, 
with two conical sections per stage. 

5.3 Applications 

The most direct application for the program is to 
refine the blade shapes of an axial turbine in preliminary 
design. However, less specific applications to research 
type problems may be fruitful. Used in conjunction with 
a sufficiently general simple radial equilibrium 
preliminary axial turbine design program, the program 
presented in this thesis could be used to generate a 
series of modified designs to explore the ramifications 
of flare and actuator disk effects. Perhaps it is 
possible to find a simple radial equilibrium tangential 
velocity distribution which will be less effected by the 
modification, or will "compensate" for the effects. 



88 



REFERENCES 



1. Carmichael, A. Douglas, Gai TuAbZne Eng-inee/Ung 

Handbook (John W. Sav/yer: Editor), Chapter 3, 

Gas Turbine Publications, Stamford, Connecticut, 

1966. 

2. Carmichael, A.D. and Pai, B.R., fhUd Vloo) d.n H-lghZij 

EZoAO-d AxZaZ-EZoto TuabZn&6, Imperial College of Science 
and Technology, Department of Mechanical Engineering, 
London, June 1968. 

3. Horlock, J.H., AxZaZ FZow Comp/i2^6oO^ , Butterworths 

Scientific Publications, London, 1958, Chapter 1. 

4. Carmichael, A. Douglas, MIT Subject 13.26 Notes, 

unpublished, 1974. 

5. Cravalho, Ernest G. and Smith, Joseph L. , Jr., 

TkeAmodijnaivZc^i, {An hUAoductionj , Holt, Rinehart and 
Winston, Inc., 1971, Chapter 11. 

6. Korlock, JbZd. , Chapter 5. 

7. IBM SijUmlUO and Si{6tm/370 FORTRAN IV Language, IBM, 

1974. 

8. Cress, Paul, FORTRAN JV li’ZZh WATFOR And (I'ATFIl^, Prentice- 

Hall, 1970. 

9. MIT Subject 13.26 Design Project, Spring 1974. 



89 



APPENDIX A 



Variable List 

The following is an alphabetical list of variables 
which appear in the computer main program and subprograms 

The units are the same as those listed in the 
Notation section. The same variable names are used to 
represent identical quantities in both the main and sub- 
programs. In addition, the subscripts are used in 
consistent manner throughout. An intermediate variable 
is one used locally in the program to simplify arithmetic 
expressions. Zero is v;ritten 0 . Unless otherwise 
indicated, variable naming follows the implicit REAL/ 



INTEGER typing 


convention. 


AREA 


annulus area between two adjacent 
streamlines 


Aj2l(I) 


constant term in the expression for 
tangential velocity 


Aid) 


coefficient of r in the expression 
for tangential velocity 


A2 (I) 


coefficient of r^ in the expression 
for tangential velocity 


B(I) 


coefficient of 1/r in the expression 
for tangential velocity 


Cl 


intermediate gas property constant 
used by RNEXT 


C2 


intermediate variable 


C3 


intermediate variable 
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D 


intermediate variable 


DAL P HA 


stator tuning angle 


DBETA 


rotor turning angle 


DD (I, J) 


actuator disk correction factor 


DISK(I) 


function subprogram which calculates 
t'he non-dimensional actuator disk 
correction factor 


DM 


mass flow per channel between 
streamlines 


DR 


difference between the calculated 
and actual tip radius 


DROLD 


value of DR in last iteration 


DRZ 


difference between R(I,J+1) and RGUESS 
calculated by STREAM 


DRZOLD 


value of DRZ in last iteration 


F(I,1) 


total absolute gas velocity 


F(I,2) 


total relative gas velocity 


F(I,3) 


alpha; absolute flow angle 


F(I,4) 


beta; relative flow angle 


FF (I, J) 


flare correction factor 


FLARE ( I ) 


function subprogram which calculates 
the non-dimensional flare correction 
factor 


FLOW 


total turbine mass flov-7 


FVT(I) 


function subprogram which calculates 
tangential velocity as a function of 
radius 


FVX(I) 


function subprogram which calculates 
axial velocity as a simple radial 
equilibrium function of radius 


G 


g = 32.1740 f t-lbm/lbf/sec ^ 


GAMMA 


Y = ratio of specific heats; input 
constant 
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HR 


mean rotor height 


HS 


mean stator height 


I 


subscript which identifies a transverse 
station along the turbine axis 

1=1 stator leading edge 
1=2 rotor leading edge 
1=3 rotor trailing edge 


ID 


user determined input data identification 
number 


IREAD 


input device channel number 


I WRITE 


output device channel num.ber 


J 


streamline index; J = 1,2,3, ..., JJ, 
starting at the hub 


JJ 


quantity of streamlines considered 


JMl 


quantity of flov; channels considered 


K 


key variable; K = 1 modified simple 

radial equili- 
brium data 

K = 2 SRE data 


KK 


key variable; KK 1 program diagnostics 

not printed 

KK = 1 diagnostics 
printed 


LI 


loop count; number of times all the 
strearalines for the turbine are 
calculated 


L2 


loop count; number of times a set of 
streamlines must be calculated at a 
station to satisfy the geometry 


L3 


loop count; the total number of times, 
during the calculation of a set of 
streamlines at a plane, that the 
subprogram RNEXT is called 


MM 


key used in labeling output data 


MV 


absolute mach number (REAL variable) 
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MW 


relative mach number (REAL variable) 


PHIHUB 


hub flare angle 


PHITIP 


tip flare angle 


N1 


test variable for Ll loop divergence 
check 


N2 


test variable for L2 loop divergence 
check 


N3 


test variable for L3 loop divergence 
check 


PI 


3.14159 


P0 (I) 


stagnation pressure 


R(I, J) 


radial location of a streamline 


REF 


reference radius for flare correction 


RGAS 


specific gas constant; input data 


RGUESS 


temporary storage for a trial R(I,J+1) 


RHO 


static density 


RNEXT ( I ) 


function subprogram which calculates 
a trial value of R(I,J+1) 


RPM 


turbine rpm; input data 


RR 


intermediate variable used as the 
argument for FLARE (I), FVX(I), FVT(I) 


RT(I) 


tip radius 


Rj3 (I) 


input reference radius at which VX(I,J) 
VX^(I) 


S 


sum which accumulates the difference 
between streamline positions on the 
current iteration (Ll) and the 
previous (Ll-1) 


SOLD 


value of S from previous iteration 


SI 


term in S 


TESTl 


test value to terminate loop Ll 
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TEST2 


test value to terminate loop L2 


TESTS 


test value to terminate loop L3 


TEST4 


test value to ensure the axial 
velocity is positive 


TGNT 


tangent of the flare angle 0 


TR 


ratio of static to total temperature 


T0 (I) 


stagnation temperature 


Tl 


input percent used to calculate TESTl 


T2 


input percent used to calculate TEST2 


T3 


input percent used to calculate TESTS 


U 


blade speed 


VRATIO 


Vi/Vj 


VT(I, J) 


tangential velocity 


VTZ 


average value of 'tangential velocity 
between two streamlines 


VX(I, J) 


axial velocity 


VXR(I) 


reference value of axial velocity; 
input data 


VXSRE (I, J) 


simple radial equilibrium axial velocity 


VXZ 


average value of axial velocity between 
two streamlines 


VX0 (I) 


input reference value of axial velocity 


WE 


entrance width dimension used to find 
the flare correction factor 


WR 


rotor width 


WRATIO 




WS 


stator v/idth 


Z 


intermediate variable 


Zl, . . . , Z6 


intermediate variables 
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APPENDIX B 



Progx'am Listing 

This appendix contains the computer program listing 
of the MAIN program, subroutine STREAM, and function 
subprograms RNEXT, FLARE, DISK, FVX, and FVT. Refer to 
the table belov/. The purpose of each program segment 
is indicated; flow charts for MAIN, STREAI^l, and RNEXT 
appear in Chapter 3. 



paragraph 


segment 


page 


B.l 


MAIN 


95 


B. 2 


STREAM 


103 


B. 3 


functions 


106 



B.l MAIN 

The MAIN program reads the input data, governs the 
iteration loops to generate the streamlines, calculates 
the velocity triangles, and prints out the results. The 
listing begins on the next page. 
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B.2 Subroutine STREAf-1 

STREAI'i calculates an j+1' iterating on assumed 

values for V and Vq until r. and the average 

X 0 1/3+1 

properties based on r . . and r. . , are consistent. 
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B.3 Function Subprograms 



The listing for each of the five function subpro- 
grams appears on one of the following pages. 



Function 

RNEXT 



FLARE 



DISK 



FVT 



FVX 



Page Description 

107 given a value of r. . and the 

1 ,: 

properties V^, Vg, RNEXT 

calculates a trial r. 

1 , 1+1 

108 given r^ j+1' FLARE calculates a 

non-dimensional factor which, 
multiplied by V in STREAi-l, 
gives a modified for FLARE 

109 given appropriate values of 

V . . , DISK calculates a non- 

XI, 1 

dimensional factor used to 
modify axial velocity for 
actuator disk effects 

110 given r^ j, FVT calculates a Vg 

by a proposed simple quadratic 
relation 

111 given r. ., FVX calculates a V 

using a simple radial equili- 
brium relation 
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APPENDIX C 



Specific Instructions For Use of the Program 
C.l General 

This appendix is intended to provide detailed 
instructions for the use of the computer program 
developed in this thesis. Preparation of the input 
deck and interpretation of the output will be covered. 

The program was run on the Interdata 70 computer 
at the Joint Mechanical and Civil Engineering Computer 
Facility at MIT, but should be easily adaptable to any 
FORTRAN IV capable system. References 7 and 8 served 
as the FORTRAN IV language guide. All READ and WRITE 
instructions are written with dummy variables IREAD and 
IWRITE, respectively. Therefore, the only tw’o FORTRAN 
statements which need to be changed in order to run on 
another system using different input/output' channel 
numbers are the two assigning integer values to IREAD 
and IWRITE. The control cards, listed below, which may 
be system dependent, are marked by an asterisk. 

C.2 Loading Sequence 

Listed belov/, in the sequence that they should be 
loaded, are the contituents of the program: 

* (1) Job card 

*(2) Language and other compilation control 
cards 
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(3) MAIN program 

(4) Subprograms (in any order) 

a) subroutine STREAM 

b) function subprogram RNEXT 

c) function subprogram FLARE 

d) function subprogram DISK 

e) function subprogram FVT 

f) function subprogram FVX 
*(5) Execute card 

*(6) Common length designation and the 
execution control cards 

(7) Data cards 

*(8) Job termination card 

C.3 Job Control Cards* 

The program segments marked by an asterisk may 
vary in content and format from system to system or with 
time, and must be determined from an appropriate user's 
manual. It should be noted that the amount of blank 
common utilized by the program depends upon the dimension 
of some of the variables, which relates to the number of 
streamlines considered. 

C.4 Data Cards 

The data deck is composed of any number of separate 
data sets stacked directly together. Each set is 
independent of every other set and represents a separate 
simple radial equilibrium turbine design to be processed 
by the program. 
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Each data set consists of ten separate cards v/hich 
must be present and in the proper sequence. A specific 
format is given for the data on each card, although free 
format (data fields separated by commas) may be used if 
the computer system has that feature. 

The context of each card is given belov; in 
paragraph C.4.1. A superscript refers to an explanatory 
note in paragraph C.4.2. 

C.4.1 Format Table 



Columns 


Format 


Information 




Units 


CARD ONE 










1-4 


14 


number of streamlines^ 




none 


5-8 


14 


2 

diagnostics print key 




none 


9-12 


14 


user data set number^ 




none 


CARD TWO 










1-10 


FIO. 4 


gamma (ratio of specific 
heats) 




none 


11-20 


FIO. 4 


specific gas constant 


ft 


lbf/lbm/°R 


CARD THREE 










1-9 


F9.3 


turbine rpm 




rpm 


10-18 


F9. 3 


total mass flow rate 


Ibra/sec 


19-27 


F9.3 


total temperature, plane 


1 


°R 


28-36 


F9. 3 


total pressure, plane 1 




psi 


37-45 


F9. 3 


total temperature, plane 


2 


°R 


46-54 


F9.3 


total pressure, plane 2 




psi 


55-63 


F9. 3 


total temperature, plane 


3 


°R 


64-72 


F9.3 


total pressure, plane 3 




psi 
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CARD FOUR 



1-10 


FIO, 6 


hub radius, plane 1 


ft 


11-20 


FIO. 6 


tip radius, plane 1 


ft 


21-30 


FIO. 6 


hub radius, plane 3 


ft 


31-40 


FIO. 6 


tip radius, plane 3 


ft 


41-50 


FIO. 6 


stator width 


ft 


51-60 


FIO. 6 


rotor width 


ft 


CARD FIVE 








1-10 


FIO. 6 


4 

test percent 1 


none 


11-20 


FIO. 6 


test percent 2 


none 


21-30 


FIO. 6 


test percent 3 


none 


31-40 


FIO. 6 


test value 4 


ft 


CARD SIX 








1-10 


FIO. 6 


SRE reference radius, 
plane 1 


ft 


11-20 


FIO. 6 


SRE reference velocity, 
plane 1 


ft/sec 


21-30 


FIO. 6 


SRE reference radius, 
plane 2 


ft 


31-40 


FIO. 6 


SRE reference velocity, 
plane 2 


ft/sec 


41-50 


FIO. 6 


SRE reference radius, 
plane 3 


ft 


51-60 


FIO. 6 


SRE reference velocity, 
plane 3 


ft/sec 


CARD SEVEN 








1-10 


FIO. 5 


5 

B for plane 1 


f t^/sec 


11-20 


FIO. 5 


A^ for plane 1 


ft/ sec 


21-30 


FIO. 5 


for plane 1 


-1 

sec 


31-40 


FIO. 5 


A^ for plane 1 


ft sec 
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CARD EIGHT 



(Same as card seven, but for plane 2) 

CARD NINE 

(Same as card seven, but for plane 3) 

CARD TEN 

(The tenth card is blank) 

C.4.2 Notes 

1. The limit on the number of streamlines is only 
a practical one, related to machine computing speed and 
storage, and can be any reasonable number. The program 
was originally run with eleven; any greater number v;ill 
require a nev; DIMENSION statement for the main and all 
subprograms. Additional blank common space may also 
have to be allocated. 

2. If this data bit=l, all available diagnostics 
will be printed out. See paragraph C.5.3, below. If not 
equal to unity, no diagnostics except error messages will 
be printed out. 

3. The user may use any three digit number to 
identify particular data sets. Some num.ber must be 
supplied . 

4. Refer to paragraph 3.3.13. 

5. B, A^, A^, are the coefficients of powers 

of r in the tangential velocity relation. See 
paragraph 3.3.6. 
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C.5 Output Discription 
C.5.1 General 

The output for each data set may be considered in 
three parts: 

(1) introductory 

(2) diagnostic (optional) 

(3) velocity triangle data 

(i) modified simple radial equilibrium 
values 

(ii) SRE values. 

C.5. 2 Introductory Output 

Input data and other selected calculated values, 
which do not change during the execution of a data set, 
are printed out in the introductory block. The symbols 
appearing have the same meaning as used in the program 
and given in Appendix A. If the data is input also, 
the output units are the same. TESTl, TEST2, TESTS, 

REF, and WE have units of feet, and PHITIP and PHIHUB 
are in degrees. 

C.5. 3 Diagnostic Output 

If input variable KK=1, certain lines of information 
are printed out during execution of the program iterations. 
Sample output of this type is shown at the end of this 
appendix. Each line contains the current values of 
certain important variables of the program at a certain 
step in the calculations. This information was used in 
the development of the program and should not be 
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required routinely. However, it can increase confidence 
in the results, and may be useful if unexpected problems 
occur during operation in the future. Additional details 
on the meaning of the loop counters can be learned from 
the flow charts presented in Chapter 3 . 

C.5.4 Velocity Triangle Output 

A group of results is printed for each of two cases: 

(1) The streamline numbers and locations, and 
velocity results are printed for a simple radial equili- 
brium preliminary design which has been modified by the 
flare and actuator disk corrections determined by the 
program. Results pertaining to each of the planes is 
printed, followed by a block of results involving values 
from more than one plane. 

(2) For comparison purposes, values calculated 
at the same streamlines as above, but with no modifica- 
tion for flare and actuator disk are printed in the 

same format. 

The meanings of the variables are the same as 
listed in Appendix A, or are self-explanatory. The 1 
and 2 appearing in some column labels refer to rotor 
inlet and outlet triangle, respectively, following the 
usual convention. Distances are in feet, velocities in 
ft/sec, and angles in degrees. The value VXR(I) printed 
for each plane is the final value for the constant term 
in the FVX relation for axial velocity. It was changed 
from VXff(I) to satisfy continuity. 
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C.5.5 Error Messages 



The input data to this program is assumed to be 
from a feasible preliminary design. Therefore/ many 
types of execution errors, such as from the appearance 
of negative radii, should not occur. If some inconsis- 
tant or infeasible data should be used, however, measures 
have been taken to minimize the chances of hanging the 
program up in an iteration loop. The details of the 
testing procedure to accomplish this are found in 
paragraph 3.4. Each iterative loop is provided with a 
means of exit if the error term increases for at least 
two successive iterations. Another check is made to 
determine if the axial velocity is less than 1.0 ft/sec 
(an arbitrary, small positive value, to avoid the 
possibility of dividing by zero) . 

If any of the conditions listed in the next para- 
graph are discovered, an error message will be printed 
out regardless of whether or not the user has specified 
the printing of the diagnostic data. 

The follov/ing is a list of error conditions tested 
for, and their corresponding output messages; 

(1) Condition: the axial velocity computed by 

the function subprogram FVX is less than 
1.0 ft/sec 

Error message: Ll = xxx L2 = xxx L3 = xxx 

I = xxx J = xxx VX(I,J) LESS 
THAN 1.0 FT/SEC 
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ActJ-on: execution continues 

(2) Condition: the change, S, in the 1=3 plane 

streamline positions is not 
decreasing on successive iterations 
Error message: ITERATION TERMINATED DUE TO 

INCREASING VALUES OF S 
Action; execution of current data set is 
terminated 

(3) Condition: the annulus geometry error, DR, 

is not decreasing on successive 
iterations 

Error message; ITERATION TERIilNATED DUE TO 

INCREASING VALUES OF DR 
I = XXX L2 = XXX L3 = XXX 
Action: execution of current data set is 

terminated 

(4) Condition: the change, DRZ, in the calculation 

of R(I,J+1) by subroutine STREAM 
is not decreasing on successive 
iterations 

Error message: ITERATION TERMINATED DUE TO 

DIVERGING POSITION OF STREAM- 
LINE J+1 I = XXX J = XXX 

Ll = XXX L2 = XXX L3 = xxx 
Action: execution of current data set is 

terminated 



120 




f 



C.6 Additional Options 



In some applications it may be desired to find the 
simple radial equilibrium streamlines, or to consider 
the effects of the flare or actuator disk correction 
alone. To accomplish this, the executable statements 
in one or both of the subprograms FLARE (I) and DISK (I) 
may be replaced by the statement FLARE = 1.0 or 
DISK = 1.0, as appropriate. 
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Figure C-1 
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APPENDIX D 



Integration of the Simple Radial Equilibrium Equation 

A convenient form of the simple radial equilibrium 
equation is 



dV 



-V 



X 



V 



dr 



e d , ,, , 
r dr ^ 6 



(D.l) 



Given explicit functions, Vg(r), it can be easily 
integrated to find a relation for axial velocity as a 
function of radius . 

Suppose the expression for the tangential velocity, 
Vq (r) , can be written in the following form: 

VQ(r) = a2r^ + a^r + a^ + |. (D.2) 

The bottom sign applies after a stator row, and the top 
sign applies after a rotor rov;. 

Then, preceding to evaluate the right hand side of 

equation (D.l) , 

rVg = ; b 

A(rVg) = + 2a^r + a^ 

i A(rVg) = 3a2r + 2a^ b (D.3) 
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Finally, multiplying (D.2) by (D.3), (D.l) becomes 



•\-dF ' Sa^a^r^ + (4a_^a2 + 2a^Mr 

, a b 

+ (3a^a. + 3a_b) + (a ^ + 2a.b)i + -2- (d.4) 

01.2 o 1 r ^2 

Equation (D.4) may be integrated directly from a 
reference radius, r^, and axial velocity, to an 

arbitrary radius r, and corresponding V . This gives 

- - ^xi') = - ^i') - r.^) 

+ ( 2 a^a 2 + a^^) (r^ - r^^) + 3(a^a^^ + a^b) (r - r^) 

+ (a ^2 - 2a^b)ln(|-) ± a^b (^ - ^) . 

Rearranging, 

- (4a^a2 + 2 a^M(r 2 - r^M _ 6 (a^a^ + a 2 b) (r - r^) 

- (2a ^ + 4a b)ln(r/r.) + 2a b(p - ^) • (D.5) 

The general expression, (D.5), may be simplified 
for two special cases in common use: 
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(i) Let 32 = = 0; , b 0 . Then (D.5) becomes 

^ _2a^Mr' - ± 4a^b ln(r/r^). 

(D.6) 

This represents the constant reaction type design.^ 

(ii) Let a„ = a^ = 0 ; a / b 7 ^ 0 . Then (D.5) becomes 



V ^ - V 

X XI 



-2a ^ ln(r/r.) + 2a b(i - ~ 
o 1 o r r 



(D.7) 



This represents a design in which the stator inlet angle, 

g 

is approximately the same at all radii. 
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APPENDIX E 



Actuator Disk Superposition 

Approximate solutions of the actuator disk theory, 
discussed in Chapter 2, are, in the case of the 
isolated disk: 



V 

X 



(V ) 

X 



(V ) 

X 






,TTX, 

exp(-pr) 



(E.l) 



upstream of the disk (x < 0) , and 



V 

X 









, irx, 
exp(— 



(E.2) 



dovmstream of the disk (x > 0) . 

These equations pertain to the axial variation of 
axial velocity along a streamline. 

Suppose it is desired to consider the effects of 
these actuator disk solutions at the three principal 
planes of a single stage axial turbine. Two actuator 
disks are involved, corresponding to the stator row and 
the rotor row. The three planes of interest are the 
stator leading edge (plane 1), the rotor leading edge 
(plane 2), and the rotor trailing edge (plane 3). 

For simplicity, assume that the stator leading 
edge is influenced only by the stator actuator disk 
(for equal blade widths, the rotor disk is three times 
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the distance away) , and that the rotor trailing edge is 
influenced only by the rotor actuator disk. The rotor 
leading edge feels the effects of both the stator and 
rotor disks. Refer to Figure E-1. 

The relations (E.l) and (E.2) apply to the stage 
inlet and outlet planes respectively. The effect at 
the rotor leading edge, point S, can be obtained by 
superimposing the velocity corrections due to the two 
disks. The effect of the stator upon S is given by 
(E.l) v;ith X = -WS/2 and h = HS ; the effect of the 
rotor is given by (E.2) with x = +WR/2 and h = HR. 

In addition, the velocities at infinity, appearing 
in the two equations must be interpreted. Let (V ) 

I 

denote the axial velocity a long way downstream of the 
rotor disk; similarly, let (V^^)_^ denote the axial 
velocity a long way upstream from the stator, and so 
on. Then (E.l) and (E.2) become 






(E.3) 



2 



which is the effect of the rotor, and 





(V ) 
xs 






, Tr‘WS 



(E.4) 



2 



which is the effect of the stator upon point S. 
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Figure E-1 



128 



Now assume 



^ simple radial equili- 

brium velocity at the 
rotor leading edge 
(plane 2) corrected for 
flare; 



(V ) _^ = V = the axial velocity at plane 1; 

X S X ““ 



(V ) . oo “ ^ 4 . ~ axial velocity at plane 3. 

X XT • X i 



Then finally, the corrected axial velocity, '^^ 2 ' 
point S (plane 2) can be v;ritten 



V 

V 



x2 



= 1 + 



V - V V 

x+ xm , TT*WR, xm 
exp (-^— T — ) - 



- V 



xm 



2V 



xm 



2 ‘HR' 



2V 



xm 



X- , TT*WS. 



(E.5) 

V ^2 the axial velocity corrected for both flare and 
actuator disk effects. 
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APPENDIX F 



Calculation of the Flare Geometry Factors 

The pii.rpose of this appendix is to derive the two 
parameters REF and WE used by the computer program to 
compute the flare angle <}). The computation of <p at an 
arbitrary point, Q, within the annulus is also demon- 
strated. The stage inlet and outlet hub and shroud 
radii are given information. 

Referring to Figure F-1, on the next page, it is 
seen that the point P is the intersection of the extended 
line segments forming the hub and shroud walls of the 
turbine. The line x = 0 represents the inlet plane of 
the turbine, x = x^ represents the outlet plane, and 
y = 0 represents the axis of the machine. 

The equations for the lines representing the shroud 
and hub walls, respectively, are v;ritten as 



Y = 



*”t3 *’tl* ^ „ 

- — ^ 



Y = 



‘■'hs - ■'hl> 



X. 



X + R. 1 . 
hi 



(F.l) 



(F.2) 



Solving (F.l) and (F.2) simultaneously for x and y, the 
coordinates of the point P, namely (-WE,REF), are found 
to be 
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y 
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WE = 






(WE > 0) 



(F. 3) 



WF 

= ■’hi * "*hl " "h3> 



(F.4) 



The angle 4> for some arbitrary point, Q, within the 
annulus, can be found from the expression 



tan (p 



y - REF 
-‘o 

WE + X 

o 



(F. 5) 



where the coordinates of the point 0 are (x , y ) . 

o o 
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